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Abstract. We discuss the effect of the stationarity on the avalanche statistics of 
Barkhuascn noise signals. We perform experimental measurements on a FcgsBis 
amorphous ribbon and compare the avalanche distributions measured around the 
coercive field, where the signal is stationary, with those sampled through the entire 
hysteresis loop. In the first case, we recover the scaling exponents commonly observed 
in other amorphous materials (r — 1.3, a — 1.5). while in the second the exponents 
are significantly larger (t — 1.7, a — 2.2). We provide a quantitative explanation of 
the experimental results through a model for the depinning of a ferromagnetic domain 
wall. The present analysis shed light on the unusually high values for the Barkhausen 
noise exponents measured by Spasojevic et al. [Phys. Rev. E 54 2531 (1996)]. 
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1. Introduction 

Crackling noise is the typical response of many complex systems displaying avalanche 
dynamics 'V . When perturbed by a slow varying external force or field, these systems 
respond with impulsive events which span a broad range of sizes. This noise occurs in 
micro- and meso-scopic systems such as superconductors, ferromagnets, charge density 
waves, fluids in porous media, fracture, or plasticity, but also on large-scales, as in the 
tectonic motion of faults responsible of earthquakes. 

The existence of a such a wide class of phenomena displaying similar noise features 
could imply that some general aspects of the dynamics should always be present, 
irrespective of the microscopic details of the system. One possible explanation is that 
all these systems are in the proximity of some non-equilibrium critical point and for 
this reason they display power law distributions over several decades, together with 
l//"-type noise. To understand this behavior, statistical mechanics models have been 
proposed. In some case, these are able to predict the values of the critical exponents 
and provide important tools to analyze and interpret the experimental data [Tj. A 
significative example is the analysis of the average shape of the noise pulse jll |2j, 
which was predicted theoretically to be asymmetric, while experiments, for instance 
in magnetic systems, show instead a marked leftward asymmetry This discrepancy 
has been recently resolved, showing that the asymmetry is due to the negative effective 
mass of the magnetic objects (i.e. the Bloch domain walls) responsible for the noise 
|3]. This example illustrate how simple models can improve our understanding of the 
magnetization dynamics. 

A very important feature of crackling noise is the condition of stationarity during 
the dynamics. In fact, there are systems which typically respond with avalanches during 
a transient regime, the clearest example being the acoustic emission prior to fracture 
in brittle materials [Sj, but also the magnetic noise emitted during the fracture of steel 
[H]. In other cases, it is not always obvious to understand wether the system is in a 
stationary state or not. The question is important because non-stationarity can bias 
the scaling of the avalanche distributions. A notable theoretical example in this respect 
is offered by the Random Field Ising Model (RFIM) [7j, where the critical exponent of 
the avalanche size distribution r, changes in three dimensions from r = 1.6 to r = 2 in 
stationary and non-stationary conditions, respectively. The reason behind this result is 
that in non-stationary condition the distribution is integrated over different values of 
the control parameter, yielding a larger effective exponent. This mechanism was also 
discussed in general terms by Sornette in Ref. P|. 

Experimentally the problem of stationarity has been addressed with a certain detail 
in the context of Barkhausen noise, the crackling noise occurring in ferromagnets as the 
magnetic field is increased. In 1981, Bertotti and Fiorillo P suggested to investigate 
the noise properties only in a limited range of magnetization curve around the coercive 
field, in order to ensure the stationarity of the detected noise signal. In this regime, 
the magnetization process is ruled by the jerky motion of domain walls, the interfaces 
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delimiting two domains with opposite magnetization. Before this fact was pointed out, 
noise experiments had been performed all along the hysteresis loop, i.e. in highly non 
stationary conditions, leading to results sometimes difficult to interpret theoretically. 
This new approach to the investigation of Barkhausen noise generated a series of 
important studies which have highly improved our undestanding of the phenomenon 
jlOj . Nevertheless, this scientific field is still very active, also because the Barkhausen 
effect is generally considered as one of the cleanest examples of crackling noise pp. 
A detailed analysis of the vast scientific literature on the subject can be found in our 
recent review ^T]. The main conclusion is that, despite the large variety of soft magnetic 
materials, most of the experimentally measured scaling exponents can be grouped into 
two different universality classes, reflecting the dominating magnetic interactions |12j . 
Nano crystalline and polycrystalline materials are dominated by long range interactions, 
originated by stray flelds inside the sample. Amorphous materials, especially when 
subjected to an external tensile stress [12], are dominated by the short range elastic 
tension of the domain wall. 

As a matter of fact, there is a number of signiflcative experimental results which do 
not flt into these two classes. A few materials show critical exponents which are within 
the values found for the two classes, especially in systems where the domain patters 
is highly fragmented, and dominated by residual stresses. This is the often the case of 
unstressed amorphous alloys, where the mechanism of production leaves a high degree of 
disorder and quenched-in stress in the material. Thus one may consider these variations 
as experimental fluctuations. 

A more serious problmem comes from the data reported in the paper of Spasojevic 
et al. ^1], which also employed an amorphous sample in an unstressed state, but where 
the critical exponents are signiflcantly larger than the ones found in the two classes. 
This result has often been referred toas the signature of another universality class in 
soft magnets, although no convincing theoretical explanation has been proposed. In 
our opinion, it is difficult to believe that this new universality class exists, as it would 
require that in the sample employed in Ref. jHj domain wall dynamics would be ruled by 
completely different interactions than in other apparently similar amorphous alloys. We 
are more inclined to think that in Ref. pi] the experimental conditions are significantly 
different from the standard ones. Assuming that the noise signal was not stationary, we 
can try to explain the observed high values of critical exponents. In fact, very similar 
values of the exponents were measured earlier in NiFe by sampling the noise over the 
entire hystresis loop ^S]- On the other hand, following the recommendation of Ref. [H], 
the authors of Ref. ^1] actually measured the signal in a small magnetization window 
around the coercive field. While this is certainly true, we will show in this paper that 
non stationarity can be intrinsic to the experimental setup and material properties. 

To corroborate our hypothesis, we perform a series of Barkhausen noise 
measurements on a similar amorphous sample in non stationary conditions, obtaining 
critical exponents very close to the ones reported in Ref. Moreover, we also 

measured the exponents under stationary conditions, recovering the values typical 
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Figure 1. Hysteresis loop of an FegsBis amorplious alloy under moderate tensile stress 
(10 MPa). The signal detected between A and A', corresponding to a magnetization 
range of about 1.8 T, is highly non stationary, as revealed by the corresponding time 
signal in the bottom right corner. The signal between B and B' is stationary, and 
corresponds to a range of about 0.1 T. 

of amorphous alloys [ini 1121 HE]- In addition, we show that there exists at least a 
single critical exponent which is not affected by non-stationarity. This exponent, not 
determined in Ref. ^1], relates the average size of an avalanche to its duration, and 
was shown in |2] to coincide with the high frequency exponent of the power spectrum. 
Hence, we verify that the noise spectrum does not depend on the stationary state of the 
dynamics. Finally, we show that a model of domain wall depinning [TI) t ll7 | ni ^ nT H I ^ I^ 
can be used to explain quantitatively the scaling exponents for both stationary and non- 
stationary conditions. 

A similar comparison between stationary and non stationary signals can be very 
important not only for magnetic systems, but for any other complex systems showing 
crackling noise. The lesson we draw from our analysis is that one need to check carefully 
for the stationarity of the signal before a detailed comparison between experiments and 
theory can be attempted. 

2. Experimental conditions for stationarity and non-stationarity 

In our cited review jjllj, we have analyzed many experimental papers of the literature 
to establish a set of practical rules to be able to compare experimental data with the 
theoretical predictions. A fundamental condition is to utilize an experimental setup 
which ensures the stationarity of the noise signal. This requires to take care of the 
possible variations of the experimental conditions during the measurements and/or 
inside the materials. 

A crucial feature, in this respect, is the role of the demagnetizing field which is 
not constant in non-ellipsoidal samples, being is especially intense close to the ends of 
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the sample. A practical rule is thus to detect the signal in a region where this field 
is reasonably constant, as, for instance, using a pickup coil with limited width (a few 
mm, in practice). Another important aspect, as mentioned in the Introduction, is to 
delimit the acquired signal to a small magnetization region around the coercive field. In 
practice, this means to consider a region of the hysteresis loop having a roughly constant 
permeability. Neither of these two conditions are fulfilled in the paper of Spasojevic et al. 
In fact, their amorphous sample is pretty short (4 cm), and thus the the demagnetizing 
field has a strong spatial variation which highly affects the magnetization dynamics. 
In addition, the pickup coil is much larger than the sample length: this implies that 
the detected noise is the superposition of signals of the coils all along the samples, 
which correspond to cross sections of the sample having different permeabilities. This 
measurement is thus roughly equivalent to consider a material with a curvy hysteresis 
loop, and detect the signal all along the loop. 

This is exactly the experimental condition we have considered in order to reproduce 
results similar to Spasojevic et al obtained for a Vitrovac amorphous material. In 
fig. H we depict the typical curvy hysteresis loop and the Barkhausen noise in a 
highly disordered amorphous ribbon. We have used an FessBis high magnetostrictive 
amorphous alloy, having dimension of 28 x 1 x 0.002 cm. We have applied a moderate 
tensile external stress, as the unstressed sample show a poor signal-to-noise ratio, and 
does not permit to determine the critical exponents with sufficient accuracy. The 
non stationary signal is detected along a quite extended part of the hysteresis loop, 
corresponding to about 1.8 T (A-A' in fig. Q), and the applied stress is about 10 MPa. 
This region is visibly much larger than the linear region of constant permeability around 
the coercive field. Thus the stationary signal is detected in a much smaller magnetization 
region (B-B' in fig. Q]). To improve the signal-to- noise ratio we have slightly increased 
the tensile stress to 50 MPa. 

The details of the experimental setup have been described at length in other papers 
O HI]. Essentially, a 30 cm solenoid provides a constant applied field all along 
the sample. We apply a triangular field with peak amplitude of 70 A/m, reaching a 
magnetization of 1 T, smaller than the saturation value. The coercive field is about 6 
A/m. To resolve a well defined sequence of avalanches, the frequency of applied field is 
around 10-20 mHz. The tensile stress is applied using weights of 200 and 1000 g (about 
10 and 50 MPa, respectively). 

3. Experimental results 

The statistical properties of the Barkhausen noise are usually characterized by a set 
distribution functions, typically displaing scaling. Estimating several scaling exponents 
in a single material has highly improved the possibility to verify the theoretical 
predictions and compare different systems. In the present experiment we will try 
establish what is the set of exponents that depends on the condition of stationarity. 
The first step to calculate the distributions is to extract the sequence of avalanches 
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avalanche duration 7(s) avalanche size S (Wb) 



Figure 2. Avalanche duration and size distributions for a stationary (blue) and non- 
stationary (red) signal. Critical exponents a and r of Eas. lll2l are fitted in the linear 
part of the plots. Two values of the applied field frequency (lOmHz and 20niHz) are 
reported, showing that the distributions do not depend on the driving rate. 



(or pulses) from the noise signal, introducing a small threshold which defines the 
beginning and the end of each avalanche. The distance between these points is defined 
as the avalanche duration T, while the size s is the time integral between the same 
points. The identification of a small threshold is a procedure which works pretty well in 
materials where the signal-to-noise ratio is relatively large, as for instant in FeCo based 
amorphous alloys under tensile stress jUlEI- For more disordered materials, or with a 
maze domain structure, as the one used here, the estimation of the critical exponents 
is a bit less accurate because of the reduced signal-to-noise ratio. In any case, we are 
interested to see wether there is a significative change of the values of the exponents 
beyond their uncertainty level. 

We first consider avalanche duration and size distributions scaling as 

P(r) = T-^g{T/To) (1) 

and 

P{S) = S-^fiS/So) (2) 

where Tq and Sq are the cutoff values, and a and r the critical exponents. In Fig. |21we 
plot these distributions for both non-stationary and stationary signal, together with an 
estimate of the critical exponents obtained by fitting the linear part (in the log-log scale) 
of the data. For the non- stationary signal, we get a = 2.1 ± 0.1, and r = 1.70 ± 0.05. 
The larger error for the duration distribution is justified by the limited range of the 
data (about two decades) and the influence of the background noise in the estimation 
of shortest avalanches. These values are very close to the ones reported in ref. as 
a = 2.22 ± 0.08, and r = 1.77 ± 0.09. For completeness, we have also considered the 
energy distribution P{E) ~ -E'"'^, integrating over time of the squared signal, getting 
e = 1.48 ± 0.05, comparable with the value od 1.56 ± 0.05 obtained in ref. [T^ . 
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For the stationary case, we get a = 1.65 ± 0.08, and r = 1.38 ± 0.04. These values 
are shghtly larger than the values reported for materials belonging to the short range 
class, where a ~ 1.5, and r ~ 1.30. It is likely that this signal is only nearly stationary, 
and a small contribution of non-stationary avalanches change a little the distribution. 

We have also considered another important critical exponent, which relates the 
average avalanche size to its duration 



where a is the exponent for the avalanche characteristic size, z/ is the correlation length 
exponent and z is the dynamic exponent 0. The exponent l/auz was shown as well 
to describe with good accuracy the high frequency limit of the power spectra [2 El CI]- 
In fig. El we plot both this distribution and the power spectra for better comparison in 
case of stationary and non- stationary signal. Remarkably, the exponent l/auz is the 
same and coincides with the theoretical value of 1.77 expected for amorphous materials 
|llj . At the same time, the two power spectra yields at high frequency with the same 
critical exponent. The power spectra have been normalized by the value of average 
induced flux rate, proportional to the average magnetization change in the material. 
The differences at low frequencies are simply due to the different temporal correlations 
between avalanches, which clearly change in the case of non-stationary signal. 

Unfortunately, Spasojevic et al. did not actually estimated the exponent l/auz, 
even if they plotted the joint area-duration distribution. They observed that this 
distribution is bounded by two lines with exponents 1.3 and 1.63 (see Fig. 6 of Ref. [E]), 
but it is difficult to estimate the exponent l/auz from this. 

4. Theoretical considerations and numerical simulations 

In order to understand the behavior observed in experiments, we consider a typical 
domain wall model (TEl HZl dl UH 1201 lUj and discuss how the non-stationarity of 
the noise signal affects the resulting statistics. We consider a flexible 180° domain 
wall separating two regions with opposite magnetization directed along the z axis. In 
absence of surface overhangs, we can describe the position of the domain wall by a single 
valued function h{f,t) . The dynamics of the domain wall is determined considering 
the contribution to the local force due to exchange and magnetocrystalline interactions, 
magnetostatic and dipolar fields, and disorder The resulting overdamped equation 
of motion is given by UHl El dl 113 1201 IHl 



where F is the eddy current damping constant, H is the applied field increasing at 
constant rate, k is an effective demagnetizing factor, h is the center of mass of the 
wall, 7^ is the domain wall surface tension and r] is an uncorrelated Gaussian pinning 
field. Eq. (H} has been shown to quantitatively reproduce the statistical properties of 
the Barkhausen noise in amorphous alloys under stress [TSl 1121 • 




'1/(TI/Z 



(3) 




H -kh + -i^V^h{f, t) + +r/(f, h) 



(4) 
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Figure 3. (left) Average size of a avalanche as a function of its duration for 
stationary and non-stationary signals. From this plot, it is possible to estimate the 
exponent 1 /avz of eq. |31 We also plot the theoretical estimation using the short range 
universality class {l/avz ^ 1.77). (right) The power spectra of the two signals coincide 
at high frequency, and are well approximated by j^'^l'^^^ 



The observed scaling behavior is related to the underlying domain wall depinning 
transition, which strictly speaking only occurs for = 0. In this case the wall is pinned 
unless the external field overcomes a critical field iJc? above which a steady motion 
ensues. In general, for interfaces close to the depinning transition, the response to small 
variations of the applied field occurs by avalanches whose sizes S are distributed as 

P{S)^S-^!{S[E,-Hfl% (5) 

indicating that the characteristic avalanche size 5*0 diverges at the transition with an 
exponent a. Similarly the distribution of avalanche durations scales as 

P{T) ~ T-'^giTiE - H.y/^). (6) 

The critical behavior associated to the depinning transition has been studied using 
renormalization group methods \2'2\ EHl 123 ES] and the avalanche exponents can be 
obtained by scaling relations ^Hl- In particular the renormalization group predicts 
r = 1.24 and a = 1.51 [TT] from two-loop e = 4 — d expansion to order O(e^) [2S] while 
simulations yield r = 1.27 and a = 1.51 . The cutoff exponents are approximately given 
by 1/a = 2.2 and 1/A = 1.3 respectively JI]. Finally, the power spectrum exponent is 
7 = 1.77, equivalent to exponent ruling the scaling of the size of an avalanche with its 
duration [TT] . 

The discussion above refers to the avalanche statistics sampled at a constant 
applied field close to the depinning transition, but this rarely corresponds to normal 
experimental conditions for Barkhausen noise measurements (for a notable exception 
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see the experiments reported in Ref. [2^1 for thin films). Typically the field is ramped at 
constant rate and k > 0, yielding an effective field -f^e// = ct — kh acting on the domain 
wall. In these condition, if H^ff < He the domain wall is pinned and the effective field 
grows until -ffe// > when the domain wall start to move increasing the restoring force 
provided by the demagnetizing field. Thus the effective field remains close to depinning 
transition and the domain wall displays a stationary avalanche dynamics, with exponents 
r, a, 7 derived above and cutoffs depending now on k ^21- These conditions are met 
for the stationary signal reported here, as well as for many other experiments reported 
in the literature for amorphous alloys ^1], and indeed the exponents are in perfect 
agreement with the theory. 

If one would ramp the field slowly and sample the avalanche distribution over all 
the values of the effective field, the result will be different. This condition would occur 
if we set = and increase the field up to ifc, but also for /c > if we consider the 
transient regime before the steady state is reached. To describe this non-stationary 
regime, we need to integrate Eq. over H obtaining 



with Tint = r + cr. Using the values of r and a reported above, we obtain Tint = 1-72. 
A similar discussion can be repeated for the avalanche duration distribution, yielding 
otint = a + A = 2.25. The exponent 7 will of course not resent by the integration and 
hence should be the same both for stationary and non-stationary conditions. 

From the discussion above we can conclude that the differences between stationary 
and non-stationary avalanche signals is simply due to the integration of the scaling 
function in the latter case. As a further illustration, we have simulated an automaton 
version of Eq. and computed the avalanche distributions in the steady-state regime 
and in the initial transient. As shown in Fig. |3J the scaling exponents reproduce 
quantitatively the experimental results. In particular, we find r = 1.3 and a = 1.5 
in stationary conditions and r = 1.75 and a = 2.25 for non-stationary conditions. 
Fig. El shows as welU that 1 / auz is does not depend on the stationarity if the signal. 

5. Conclusions 

We have measured the avalanches distributions in an amorphous material with 
high quenched-in disorder, performing a comparitive analysis of stationary and non- 
stationary signals. We have shown that the size and duration distributions are highly 
affected by non-stationarity, while the scaling relation between the average size and 
duration is a robust quantity. This result is actually non surprising, since the non- 
stationarity only changes the probability of having an avalanche (and thus the overall 
statistics), but does not affect its internal structure. 

The present discussion can be extended to other classes of crackling noise where 
problems of non-stationarity may arise. From our analysis the exponent 1/avz emerges 
as the most reliable test for the universality class of crackling noise statistics, while one 
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Figure 4. The avalanche distributions from numerical simulations performed 
in stationary and non-stationary conditions: (a) Size distribution, (b) Duration 
distribution. 
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Figure 5. Average size of a avalanche as a function of its duration for stationary and 
non-stationary signals obtained from numerical simulations. 



should be aware of the possible bias in the distribution exponents. For the Barkhausen 
noise case, we have provided a quantitative explanation of the experiments by analyzing 
a model for domain wall depinning. The exponents predicted theoretically for the model 
in stationary and non- stationary conditions reproduce to a very good precision the 
experimental measurements. A similar resoning would apply to other kind of critical 
phenomena as discussed in Ref. jHj. 

Finally, we think to have clearly demonstrated that the experiment of Spasojevic et 
al. J5 does not imply the existence of a third universality class in the Barkhausen noise, 
but only that the results were probably obtaining under non-stationary conditions. A 
similar discussion applies to the earlier measurements reported in Ref. |15j, where the 
distributions were sampled along the entire hysteresis loop. 
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